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We investigate here gravitational collapse of a perfect fluid with a linear isentropic equation of state 
p = hp. A class of collapse models is given which is a family of solutions to Einstein equations and 
the final fate of collapse is analyzed in terms of the formation of black holes and naked singularities. 
The collapse evolves from a regular initial data and the positivity of energy conditions and other 
physical regularity conditions are satisfied. As we provided here an explicit class, this gives useful 
insights into the endstates of collapse with a physically reasonable and relevant equation of state 
and for the cosmic censorship hypothesis. 
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I. INTRODUCTION 

The cosmic censorship conjecture (CCC) articulated 
by Penrose [T] is fundamental to many aspects of the- 
ory and astrophysical applications of black hole physics 
today. Despite many attempts over past decades no the- 
oretical proof or even any satisfactory mathematical for- 
mulation of CCC is available as of today in the case of dy- 
namical gravitational collapse. In the mean time, many 
authors have studied mainly spherical gravitational col- 
lapse of a massive matter cloud within the framework of 
general relativity. As the nuclear fuel of a massive star 
exhausts, it loses its equilibrium and gravity becomes the 
central dominant force which lends the star to its perpet- 
ual collapse. The gravitational collapse studies then show 
that the collapse end state is either a black hole (BH) or 
a naked singularity (NS), depending on the nature of the 
initial data from which the collapse evolves, arising from 
a regular initial state to the final super dense state [2H2] . 

The works such as above analyze physically relevant 
general matter fields that include most of the known 
physical forms of matter like dust, perfect fluids, and 
such others. Wide classes of solutions to Einstein equa- 
tions are shown in these cases to exist where the collapse 
end state is a black hole or naked singularity depending 
on the nature of initial data and the collapse evolutions 
allowed by the Einstein equations subject to regularity 
and physical reasonability conditions. 

An equation of state for the matter is, however, not al- 
ways assumed in many of these investigations. As such, 
the physical characteristics of the matter field constitut- 
ing the star are described by an equation of state relating 
the density and pressures, and therefore it is important 
to know if the naked singularities would still form when 
a suitable equation of state is assumed for the collaps- 
ing cloud. But it is to be noted that such an equation 
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of state is not precisely known for a collapsing massive 
star which will be assuming super dense states of mat- 
ter closer to the final later stages of the collapse, where 
the physical region has ultra-high densities, energies and 
pressures. It is far from clear as of now what would be a 
physically realistic equation of state that would describe 
such an ultra-high density region of collapse, and also 
whether the equation of state would remain unchanged 
or it would actually evolve and change as the collapse 
develops. 

Nevertheless it is important to understand the final 
fate of a gravitational collapse when a suitable equation 
of state is assumed, and to know if the naked singularities 
still form as collapse final states. One can thus choose 
the equation of state to be linear isentropic or polytropic 
to describe the collapse of a massive star, right after its 
departure from the equilibrium configuration where grav- 
ity was balanced by the nuclear reactions taking place at 
the core region of the star. 

The gravitational collapse of a perfect fluid with a lin- 
ear equation of state is of interest from both theoretical 
as well as numerical relativity perspectives. Our moti- 
vation also comes from certain other questions such as, 
what if the value of k increases in the range < k < 1 
when a naked singularity appears as collapse final state 
for a given spacetime dimension and for a critical positive 
tangent to the singularity curve as to be discussed below. 
In such a case will the singularity sustain its nature, or a 
formation of event horizon shall precede the formation of 
singularity to preserve the CCC? We believe as answer 
to these and similar issues would be important to under- 
stand better the physical conditions for collapse and the 
role of an equation of state towards providing a suitable 
mathematical and physical formulation of CCC. 

We therefore consider here a linear isentropic equation 
of state, p = kp, < k < 1 in the study of spheri- 
cal gravitational collapse of a perfect fluid. Self-similar 
perfect fluid collapse models with a linear equation of 
state were considered through numerical simulations by 
Ori and Piran [5] and analytically by Joshi and Dwivedi 
[6] to show how black holes and naked singularities de- 
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velop as collapse final states in this scenario. Further, 
Goswami and Joshi [10] studied the case of an isentropic 
perfect fluid with a linear equation of state in four dimen- 
sional spacetime without the self-similarity assumption, 
wherein they showed that the occurrence of BH and NS 
evolving from regular initial data depends on the choice 
of rest of the free functions available. It was also proved 
that in a general TV-dimensional dust collapse, the occur- 
rence of NS can be removed when one goes to a higher 
dimensional spacetime, thus restoring CCC [11]. Dad- 
hich et al 12 studied spherically symmetric collapse of 
a fluid with a non- vanishing radial pressure in higher di- 
mensional spacetimes with an equation of state p r = kp 
using the construction of the Joshi-Dwivedi root equa- 
tion that governs the nature (BH versus NS) of the cen- 
tral singularity. Also other cases of collapse with pressure 
have been studied by various authors to understand the 
collapse final fate [13] , 

Our aim here is to give an explicit class of collapse 
models for the perfect fluid case with a linear equation 
of state, and to analyze the same for the collapse end 
states in terms of formation of black holes and naked 
singularities. We analyze here the marginally bound case 
for the sake of simplicity and transparency and the results 
show that BH and NS form for a wide range of values of 
the parameter k in the equation of state p = kp. 

The plan of work is as follows: In section II, we study 
Einstein field equations for spherically symmetric met- 
ric. The dynamics of collapse for occurrence of NS/BH 
phases due to the effect of regular initial data is studied 
in the subsection II. 1. In subsection II. 2, we obtain the 
mass function M (r, v) as a class of solutions to the con- 
straint equation that enters through dynamical equations 
for p = kp. The implications of this function on radial 
null geodesies is investigated in the next subsection II. 3. 
Apparent horizon curve is obtained in section III and it 
is investigated for increasing value of the parameter k, 
providing illustration of its behavior in subsection III.l. 
The conclusions are specified in section IV. 

II. COLLAPSE WITH A LINEAR EQUATION 
OF STATE 

The general spherically symmetric metric describing 
spacetime geometry of a collapsing cloud can be de- 
scribed in comoving coordinates (t, r, 9, <j>) by the metric, 

ds 2 = _ e 2»(t,T) dt 2 + e 2^[t,r) dr 2 + R 2 ^ ^2 (1) 

where 

dn 2 = d6 2 + sin 2 (6)dc/) 2 

is the metric on a two-sphere. We aim here at the ques- 
tion, whether the choice of an equation of state p = kp, 
k € [0, 1] contributes through the parameter k in the de- 
velopment of BH/NS phases. The stress-energy tensor 
for the type I matter fields, and in particular for perfect 
fluids, is in a diagonal form in a comoving coordinate 



system and is given by [TJ], 

n = - P , t; = if = t+= p. 

The quantities p and p are the energy density and pres- 
sure respectively. We take the matter field to satisfy 
the weak energy condition which implies p >0; p + p >0. 
The linear equation of state for a perfect fluid is, 

p{t, r) = k p(t, r) where k € [0, 1] . (2) 

The Einstein field equations for the metric are derived 
as [TO] 

F' IF 
'' H 2 H' k /,'-'// ?!) 



v'^-^Mp)]' (4) 
R'G - 2RG v' = (5) 



G-H = l--. (6) 

Here F — F(t, r) is an arbitrary function, and has an 
interpretation of the mass function for the cloud. It gives 
the total mass in a shell of comoving radius r on any 
spacelike hypersurface t = const. The energy conditions 
impose the restriction on F, namely that F(t, r) > 0. 
In order to preserve the regularity at the initial epoch, 
we have F(ti,0) = 0, that is, the mass function has to 
vanish at the center of the cloud. Since we are considering 
collapse, we have R < 0, i.e. the physical radius R of the 
collapsing cloud keeps decreasing in time and ultimately 
it reaches R = 0, which denotes the spacetime singularity 
where all the matter shells collapse to a zero physical 
radius. The functions G and H are defined as G(t, r) — 
e-^'R' 2 and H(t,r) = e~ 2t ' R 2 . 

Let us define a new function v(t,r) by v(t,r) = R/r, 
and now we use the scaling independence of the comoving 
coordinate r to write 

R(t,r) =rv(t,r) (7) 

This gives v(ti,r) — l,v(t s (r),r) = and the collapse 
condition is now written as v < 0. The time t = t s (r) 
corresponds to the shell-focusing singularity at R = 0, 
where all the matter shells collapse to a vanishing physi- 
cal radius. The introduction of the parameter v as above 
allows us to distinguish the genuine spacetime singular- 
ity from the regular center at r = 0, with v = 1 at the 
initial epoch including the center r = 0, and v then de- 
creases monotonically with time as collapse progresses to 
the value v = at the singularity R = 0. At the regular 
center the mass function F(t, r) behaves suitably so that 
the density remains finite and regular there at all times 
till the occurrence of singular epoch. From the Einstein 



3 



equation for density we see that the mass function for 
the cloud can be written in general as, 



F(t,r) = r s M(r,v) 



(8) 



where M(r, v) is regular and continuously twice differen- 
tiable. Using equation Q in equations ([3]), we obtain 



3M + r[M >r + M ;V v'] _ M^ v 
v 2 (v + rv') k v 2 



(9) 



Then as v — > , p — > oo and p — > oo, i.e. both the density 
and pressure blow up at the singularity. We rearrange 
equation (|9| as follows, 

krM } r + [(k + l)rv' + v]M v = -3kM (10) 



The general solution of equation ( 10 ) represents many 



classes of solutions but only those classes are to be con- 
sidered which satisfy the energy conditions, which are 
regular and those which give p — > oo as v —> 0. This 
means that the energy conditions and equation of state 
p = kp isolate the class of functions M(r, v) so that 
the mass function F(t, r), the metric function v(t, r) and 
other concerned functions of v evolve as the collapse be- 
gins according to the Type I field equations. 

To obtain such a class of solutions, we consider the 
condition, 



rv 
v 



(11) 



due to which the equation ( 10 ) has a general solution of 
the form, 



M(r,v)= f(re- z ^)e- 3Z(v) 
k 



r 

where Z(v)= I — 



P + i)g(v) 



l 



(12) 
(13) 



and / is an arbitrary function. Therefore, we can choose 
a physically relevant M as given by, 



M(r,v) = M sech(re- z ^) e - 3Z{v) 



(14) 



as a solution of equation ( 10 ), where M G is a positive con- 



stant and the choice of the arbitrary function g(v) should 
be such that integrand of equation (131 should be at least 



piecewise continuous on [0, 1]. Therefore, we can write 
M (v) = M(0,v) = n e- 3Z , Mx(v) = M ' r [^ v) = 



and M 2 (v) 



M, rr (0,v) 

2! 



1! 

-5Z 



Here M\(v) — is in accordance with the requirement 
that the energy density has no cusps at the center and 
that its first derivative vanishes at the center. The den- 
sity profile for this class of models then takes the form, 



p(r, v) 



n o e- 3Z ^ v hech(re^ z ^) 
v 3 [(k + l)g{v) + 1] 

-Z(v 



[3-re- z{v >tanh{re- z ^)]. (15) 



We note that at the initial epoch, g(l) = 0, Z(l) = 0, 
and the regular density distribution takes the form 

Po (r) =M Q sech(r)[3- rtanh(r)]. (16) 

Such a density profile is decreasing away from the center 
r = 0, which is a physically reasonable feature for the 
collapsing matter cloud. For example, at the epoch v = 1, 
Po{r) — at (an approximate) value r = 3.014482. Hence 
we would take such a value of r as the boundary of the 
cloud where the energy density is zero, and where the 
interior is matched to a suitable exterior metric. 

At r = we have g(v) — and therefore e^ 3Z ^ — 
and the energy density at the center takes the form, 



-3fc 



P(0,«) 



3M Q 

v 3(k+l) 



therefore p satisfies the energy condition p > at all the 
epochs. The gradient of p(r, v) is obtained as 



-M Q e~ 4Z sech(re z ) 
v 3 [(k + l)g(v) + 1} 



{4tanh(re~ z ) - re~ z [2 tanh 2 (re~ z ) - 1]} . 

At the initial epoch v = 1, we get p^ r = at the center 
r = 0. The radial coordinate r has the range < r < 
Th which defines the interior of the compact collapsing 
object, and wherein the density profile is smooth and 
monotone decreasing away from the center of the cloud 
as one moves towards the boundary of the cloud. 



II. 1. Dynamics of collapse and the BH/NS phases 

Based on the above, we would now like to determine 
the final states for such a perfect fluid collapsing cloud 
in terms of the black holes and naked singularities. As 
we shall point out, the key factor that determines this 
final outcome is the geometry of the trapped surfaces. If 
the trapped surfaces developed early enough in collapse, 
well before the formation of the singularity, then a black 
hole is the collapse final state. On the other hand, when 
the trapped surface formation is delayed as the collapse 
evolves, due to the internal dynamics of the collapsing 
cloud, in that case the singularity is no longer fully cov- 
ered by an event horizon, and we can have light rays and 
material particles escaping away from an arbitrary vicin- 
ity of the same, thus giving rise to a naked singularity. 

To decide and understand this trapped surface geome- 
try, we need to determine the behavior of the singularity 
curve in the spacetime, which describes how and when 
different collapsing shells arrive at the zero physical ra- 
dius, thus forming the spacetime singularity. 

To determine this, firstly on integrating equation Q, 
we obtain the general metric function, 



v{r,t) = -I ln(p) + a(t) 



(17) 



where I = ^ 



We can set a(t) = 0. Let A(r, v) be 
a suitably differentiable function defined by A(r, v) t v — 
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v ' /R' , and on using equations ( 11 ) and (17), it takes the 
form 

A(r, v) „ = — 

v ' ' ' pR> 



le- z {AT -re- z [2 T 2 - 1]} 
v[l+g{v)] [3-re- z T] 



(18) 



\A(rv) 1 _ ^{4 r(r,l)-r[2r»(r, !)-!]} 

[3-rT(r,l)] 1 j 

where T(r,v) = t&nh(re~ z ) and T(r, 1) = tanh(r). 
Then on integrating equation (18), we obtain 



A(r,v) 



1 le- z {AT 



'[2 T 2 - 1]} 



v[l + g(v)\ [3-re~ z T] 



dv (20) 



We note that if we explicitly find A(r,v), we then have 
an exact solution of the Einstein field equations for the 
perfect fluid collapse case. As is known this is of course a 
complex problem, and in the present case, the arbitrary 
nature of the function g(v) would forbid such a possibility 
of writing down the exact metric for the full spacetime. 
However, we can still write the metric in the neighbor- 
hood of the center of the cloud as we point out below. 

Our aim below is to find the BH/NS phases for collapse 
for a given initial data in terms of the function M(r,v), 
and we write v' = vg(v)/r = w(r,v). The mass function 
M(r,v) satisfies the regularity condition M >r (0,v) = 
and we have A s» r closer to the center r = 0. The 
function A can be expressed as an expansion and it takes 
the form, 

A(r, v) = rA^v) + r 2 A 2 (v) + r 3 A 3 (v) + G(r 4 ), (21) 
with the first terms Ai(v) given by 



A, r (0,v) 

rj 



5 I 



-2Z(v) 



3 J v v[l + g(v) 



dv, 



2! 



A 3 (v) 



A irrr (0,v) _ ~5l f 1 
: 3! 



-4Z(v) 



9 J v v[l+g{v)} 



dv (22) 



Also, the use of definition of A(r, v) in equation ^ 
yields 



G(i,r) = d(r) e 



2rA 



(23) 



where d(r) is another arbitrary continuously differen- 
tiable function of r. Following comparison with the dust 
collapse models, we can write 



d(r) = 1 + r 2 b(r) 



(24) 



where the function b(r) basically characterizes the en- 
ergy distribution for the collapsing shells. For the sake 
of simplicity and transparency of consideration and for 
a better understanding of the model, we consider here 
only the 'marginally bound case', which corresponds to 
b(r) = 0. 



As we see now, having supplied an explicit function 
M(r, v) which is a solution to the first two Einstein equa- 
tions, that determines p(r,v), and therefore the metric 
function i>(r, v) is fully determined as above. The func- 
tion A(r, v) is also determined as above and is now a 
known function, and therefore G(r, v) is also determined. 
Finally, the metric function R is determined as we shall 
mention below. We can write the metric in the neighbor- 
hood of the center r = of the cloud as, 



ds A = -p- 2 l dt 



R' 2 

e 2rA 



dr 2 + R 2 {t,r)dn 2 



(25) 



In order to determine the nature of the singularity 
curve t s (r) which corresponds to the physical singular- 
ity at R = or v = 0, the field equation ^ can be 
written in the form 



R 2 = e 2v 



F 
R 



G-l 



(26) 



Using equations (|8| and (23), this can be expressed as, 

(27) 



M(r,v) 



?2rA 



where the negative sign is chosen since for collapse we 
must have v < 0. Integrating the above equation, we 
have 



i(r, v) = U + 



e~ v dv 



M(r,v) 



(28) 



In above equation, the variable r is treated as a con- 
stant. Regularity ensures that, in general, t(r, v) is at 
least C 2 near the singularity and therefore can be ex- 
panded around the center as, 



t(r, v) = i(0, v) + r X i{v) + -^Xi{v) + 0(r 3 ) 



(29) 



where 



Xi(v) 



dt 
dr 



X2(v) 



dr 2 



r=0 



and using equations (17) and (21), we have 



*(0,«)=t i + (3M o ) 



-3 IZ 



dv 



The time for other collapsing shells to arrive at the sin- 
gularity can be expressed by 



t s (r) =t(r,0)=ti + 



'dv 



M{r,v) 



e 2rA_ 1 



(30) 



Since the energy density has no cusps at the center, that 
means M\ = 0. The singularity curve then takes the 
form, 



t s {r) =t + rxi(0) + - X2 (0) + 0(r 3 ), (31) 
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MRC curve for various values of k 
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FIG. 1. The model realistic choice of M D is shown above. The 
allowed values for M have an absolute maximum (and also a 
local maximum) around k = 0.4. This maximum value of M 
is a necessary bound for the models for all A; € [0, 1). 



where t = t(0, 0) is the time at which the central shell 
becomes singular, and which is obtained as 

t(0,0) =t t + 3 l (H o y- 3/2 x 
2M 5 



3(1 -k) 3(fc + l) 



5- k 



3k 



The time taken by the central shell to reach the singular- 
ity should be positive and finite, and hence we have the 
model realistic condition (MRC), namely that, 



2M„ 



3(1 - k) 3(fc + l) 



5-k 



3fc 



> (32) 



and it must be finite for any k € [0, 1) 
Noting that v '(0,v) — 0, we have 



Xi{v) 



1 f 1 2e~ v A 2 



M(0,v) 



2A X 



3/2 



dv (33) 



and since A 2 (v) = 0, we have Xi(0) = f° r ai l h. Also, 
at r = 0, g(v) — and therefore Z — fcln(w) and e~ 3Z — 



-3k 



\ so we have 



-2k 



6(fc+l)' 



A],A 3 (v) = 



36(fc + l) 



Since xi(0) = 0, using expressions for M (v), A\{v) and 
A 3 (v), we compute X2(0) for an analysis of the nature of 
the central singularity in terms of its visibility or other- 
wise. We get, 



xa(«) = -(3m )' 



It then follows that, 



Xa(0) = - 



(3M P ) ; 
n 5 J 2 



where C\ 



31fc 2 + 136fc+25 
(5+3fc)(5-fc) 



and C2 



5/ 



,-5fe 



Mojv) 



•2Ai 



1/2 



—Kg 2 l M G 

3(fc + 1) + 9(1 - k) 



3fc 2 -20fc+7 



(7+/c)(7-3fc)(7+5fc)(7+9/c) 



\M 2 



2A 3 + 2A\] v 



-3k 



M (v) 



2A X 



ci 
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k+l 5-k 



are functions of k. 



3/2 



dv 



(34) 



800 c 2 I 2 
9(k + 1) 



(35) 



Since the term Mq(v)/v + 2Ai(v) > 0, and with the 
suitable continuity of the functions involved in the inte- 
grand of equation (34), the sign of the term X2( v ) to- 
tally depends on the behavior of the term M-ijv + 2A 3 . 
Therefore the positive or negative sign of X2(0), for var- 
ious values of k € [0, 1), depends on the behavior of the 
functions Mi(v), and A 3 (v). We also note that the pos- 
itive constant M Q appearing in the mass function F(t, r) 
is related with density by the expression 3M Q = p(0, 1). 
So, M Q in fact represents one-third of density of the star 
at the center taken on a cross section of an initial surface 
v = 1 (at the beginning of the collapse at the instant 
t = U). 

The significance of the term X2(0) towards deciding the 



visibility or otherwise of the final spacetime singularity 
is as follows. As seen by the equation for t s (r) above for 
the singularity curve, this term and its sign basically de- 
cides the behavior of the singularity curve away from the 
center. In particular, when the singularity curve is in- 
creasing, it is shown that the apparent horizon which be- 
gins also at the central singularity at r = 0, is increasing 
necessarily, and therefore it follows that there is no for- 
mation of trapped surfaces as the gravitational collapse 
of the matter cloud evolves, all the way till the singularity 
formation epoch. In such cL C0jSG, cLS has been shown ear- 
lier [5], the singularity is necessarily locally naked and a 
family of future directed null (and also timelike) geodesies 
come out which terminate in the past at the singularity. 



G 



1 •] i, (0) for various values of k, b(r) = 0, Mo = 1 




FIG. 2. For a fixed value of the central density corresponding 
to M = 1, the quantity x^(0) takes positive values for k — 
and also for other small values of k. Then it crosses the fc-axis 
and takes negative values as k — > 1. 



tion. On the other hand, for a black hole the trapped 
surfaces form before the singularity. Thus for a black 
hole to form we require, 

t&h( r ) <t for r > , near r = . (36) 

where to is the epoch at which the central shell hits the 
singularity. In the general case this condition is violated 
when X2 > as seen from the equation for t s (r) above. 
The apparent horizon curve that initiates at the singular- 
ity r — at the epoch to, then increases with increasing 
r, moving to the future and we have t a h > to for r > 
near the center. The behavior of the outgoing families of 
null geodesies has been analyzed in detail in these cases, 
and it is known that the geodesies terminate at the sin- 
gularity in the past, which results in a naked singularity. 
In such cases then the extreme strong gravity regions can 
communicate with external observers. 

As these aspects are already discussed earlier as cited 
above, we therefore mainly focus here on determining the 
behavior of the quantity X2 for the perfect fluid collapse 
with an equation of state p — kp, and we examine the 
impact of changing the values of the parameter k on the 
visibility or otherwise of the singularity. 



& (0) for various values of k 




FIG. 3. The nature of X2(0) is illustrated through the above 
graphs. We note that if the model begins with a higher value 
of M , higher values of k would satisfy the condition X2(0) > 0. 
Thus, it is the initial central density that can specify and 
would have ramification on the nature of singularity in terms 
of its visibility or otherwise. 



Such a singularity can also be globally visible depending 
on the nature of the mass function away from the center 
[5J. On the other hand, when the singularity curve is con- 
stant (\2 and other higher order terms are all vanishing) , 
or would be decreasing, then a black hole will necessarily 
form as the collapse final state. 

Basically, a naked singularity could occur when a co- 
moving observer (at a fixed radius r) does not encounter 
any trapped surfaces until the time of singularity forma- 



II. 2. Dust collapse 

We note that for the well-studied case of the dust col- 
lapse, we have p = at k = and the MRC reduces to 
the form M Q > 0. We have 

X 1 (0) = 0andx^ = °(0) = ^ r7 ^>0. (37) 

3M /2 

In this case a naked singularity exists for the given initial 
data of mass function. 



II. 3. The Radiation collapse and other values of k 

The radiation collapse can be considered through the 
case fe = 1/3 for which the MRC takes the form M D > 5/84 
and we have Xi(0) = and 

(3) 1 / 4 r 1 ... 121 



*r /3 (o)=- 



1 2 

M 

4 ° 



1008 



M 



6864 



(38) 



Since the bracketed term has varying sign, we find that 
X2(0) < for small positive values of M Q , and for higher 
values of M Q , X2(0) changes its sign. 

In general, to understand the nature of X2(0) in the 
interval k £ [0,1), we can obtain expressions for X2(0) 
for different values of k in order to see how it changes 
with a change in the parameter values k. The graph of 
X2(0) for various specific k values is given in Fig. 3. 



III. THE APPARENT HORIZON 

The equation of apparent horizon in a spherically sym- 
metric spacetime is given by F/R = l. Use of equations 
(j8l) and ( 14 ) give us the expression 
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IV. CONCLUSIONS 




FIG. 4. An implicit graph of X2(0) specifies the regions for 
formation of naked singularity and black hole, as determined 
by the variation in parameters M D and k. 



r 2 M Q sech(re 



-Z(u)\--3Z(u) 



It is clear that as v — > 0, we must have r->0on the ap- 
parent horizon. Therefore, the only singularity that can 
eventually be visible is that at the center of the collapsing 
cloud. Certainly, it is not possible to satisfy 1 — F/R > 0, 



i _ r 2 n o sech(re- z ^)e- 3Z ^ 



< as v — > for r > 



near the singularity and away from center, that is with 
r > 0. It follows that the region surrounding the singu- 
larity cannot be timelike, and therefore any singularity 
that might eventually form near the center with r > 
must not be visible. 



The investigation of gravitational collapse with a lin- 
ear equation of state has revealed role of the parameter 
k in deciding the nature of central singularity for the 
given initial data set of the mass function M(r,v). We 
examined here the marginally bound case with b(r) = 0. 

In particular, we gave here an explicit class of mass 
functions that satisfies the Einstein equations and the 
metric obtained represents a class of solutions for an ar- 
bitrary choice of functions g(v), representing the exact 
solutions for a perfect fluid sphere collapse with an isen- 
tropic equation of state. 

The occurrence of a locally naked singularity as col- 
lapse end state is shown for these collapse models for a 
wide variety of values of A: € [0, 1), and also formation of 
black hole is seen for a range of k values for this class of 
mass functions which are an appropriate choice of M(r, v) 
that are solutions to the Einstein field equations. 

We thus see, through an explicit demonstration, that 
both black holes and naked singularities arise naturally 
as collapse final states for perfect fluid collapse with a 
linear isentropic equation of state. The class of solutions 
we gave here have several intriguing properties some of 
which we discussed here. It is thus seen that the choice of 
k £ [0, 1) influences and contributes to decide in the for- 
mation of NS/BH phases of gravitational collapse. The 
model realistic condition discussed here emphasizes the 
role of pressure in the gravitational collapse due to which 
collapse could stop momentarily, and a bounce may take 
place. In other cases the cloud would collapse continually 
without a bounce. 

We note that the present model can be extended for a 
more general case with b(r) ^ 0, and the results on the 
same shall be reported elsewhere. 
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